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Abstract 

We give a more complete calculation of 6 — > 57 decay in 2-Higgs doublet model, including 
leading log QCD corrections from rritop to in addition to corrections from Mw to m^. 
The inclusive decay rate in the first model is found to be suppressed 23% comparing with 
the calculations without the QCD running from nitop to Mw- And the enhancement up to 
5% is found in calculations of the second model. More strict restrictions to parameters of 
2-Higgs doublet model H are found. 
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1 Introduction 



The standard model(SM) has achieved a great success recent years. However, there is still a vast 
interest beyond standard model. One of the most simple extensions of the standard model is 
to add another Higgs doublet. This may be viewed as inevitable since many attempts at going 
beyond the SM leads to an enlargement of the Higgs sector (Such as Supersymmetry Model). It 
is of interest to study the 2-Higgs doublet model, hence to constrain other models. 

The decay 6 — > 57 is one of the very useful channel for study of models beyond standard 
modelQ. Recently the CLEO collaboration has observed]^ the exclusive decay B K*'j with 
a branching fraction of (4.5 ± 1.5 ± 0.9) x 10~^. A new upper limit on the inclusive 6 — 57 
process is also obtained as B{b — » 57) < 5.4 x 10"^ at 95% C.L.0. This leads to a number 
of papers [Q regarding this decay recently. It has been argued that the experimental result is 
more close to standard model predictions and provides more information about restrictions on 
models beyond SM. This deeply depends on more precise calculations of this decay. The decay 
of 6 57 and its large leading log QCD corrections have already been calculated in many 
papers]^, |^, 0, §, 0. Due to many reasons, there are discrepancies among these papers, although 
there is not much differences in numerical results. Until recently, Ciuchini, Franco, Reina and 
Silvestrini[0] performed this calculation in three regularization schemes, to solve this problem 



completely. And also some efforts are made to give a next to leading log calculations [|TT], |T2| 
which is estimated within 20% contribution. All these efforts make it easy for calculations of 
this decay in models beyond SM. 

The radiative b quark decay in 2-Higgs doublet model has also been calculated in several 
papers]^, |13|, |T5|. It is found to be strongly QCD-enhanced. In other words, the strong 



interaction plays an important role in this decay. However, there are still some uncertainties in 
these papers. All these papers do not include the QCD running from m^p to M^y. Since the 
top quark is found to be 2-times heavier than W boson( rritop = 174 ± I0^\l GeV [|16| ), and the 
charged Higgs is also expected to be heavier (in supersymmetry model m| = M^ + m\), it needs 
a detailed calculation of this effect. 

In our present paper, by using effective field theory formalism, we recalculate the b 
sj decay in 2-Higgs doublet model including QCD running from rritop to M^, in addition to 
corrections from to rrib, so as to give a complete leading log results. In the next section, we 
first integrate out the top quark and charged Higgs, generating an effective five-quark theory. 
By using the renormalization group equation, we run the effective field theory down to the W- 
scale. In section 3, the weak bosons are removed. Then we continue running the effective field 
theory down to b-quark scale to include QCD corrections from M^y to rrif,. In section 4, the rate 
of radiative b decay is obtained. Restrictions of model parameters from experiments of CLEO 
collaboration are also given. Section 5 is a short summary. 



2 QCD Corrections from /i = rritop to fi = Mw Scale 

With two Higgs doublets, one has to avoid tree-level flavor changing neutral currents (FCNCs) 
due to neutral Higgs bosons. This is usually achieved in two ways, each involving some discrete 
symmetry, the first way (model I) is to allow only one Higgs doublet to couple to both types 
of quarks|T^. The second way (model II) is to couple only one Higgs doublet to u-type quarks 
while the other couples only to d-type quarks|T^. It is of interest to note that model II occurs 
as a natural feature in theories with supersymmetry or a Peccei-Quinn type of symmetry. 



The Lagrangian of the 2-Higgs doublet Model is 
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Where V represents the 3x3 unitary Kobayashi-Maskawa matrix, Mu and denote the 
diagonalized quark mass matrices, the subscript L and R denote left-handed and right-handed 
quarks, respectively. For model I, ^ = V2/V1; while for Model II, ^ = —Vi/v2. And Vi, V2 are 
the magnitude of the vacuum expectation values of two Higgs doublets, respectively. To keep 
explicit gauge invariance, we work in a background field gauge [p!U|. 



At first, we integrate out the top quark and charged Higgs, generating an effective five quark 

theory, introducing dimension-6 effective operators as to include effects of the absent top quark 

and charged Higgs. Higher dimension operators are suppressed by factor of p'^/mf, where 

characterizing the interesting external momentum of b quark ~ m^. For leading order of 

ml/ dimension-6 operators are good enough to make a complete basis of operators: 
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Where slDix^J)l stands for i^LDi/yJ)L + {DijSLj^vhL) and the covariant derivative is defined as 

with 5(3 denoting the QCD couphng constant. The tensor appearing in P^'^ assumes the 
following Lorentz structure, the index A ranging from 1 to 4: 

Pava ~ 9livla-i Pnva ~ 9 fieri vi , ^ 

(3) 

Pjxva 9val\xi P^i/a "^^nvarl lb- 

Then we can write down our effective Hamiltonian as 

Heff = 2^/2GFVtiyi Cr{li)0,{li). (4) 

i 

The coefficients Ci{ii) can be calculated from matching diagrams displayed in Fig.l and Fig.2. 
Calculating left hand side of Fig.l, keeping only leading orders of p^/mf , we get coefficients of 



right hand side operators: 





= Cr2^ = 




Cri 


= 0, 








Cwl ^ ^' 









(5) 

In Fig. 2, coefficients of right hand side operators are all from the finite part integrations of 
left side electroweak loops. Terms like log(/x^/mj) vanish here, because of the matching scale 
II — rrit. They will be regenerated by renormalization group running of electroweak later. After 
calculation one has 

^ + 5 ^ ^\ 1 + x x,\ 
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Cp3 = 
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Qfc Vt;iy V2(l-a;)3 (1 - x)^ 
Where 5 = M^/rn^, x = rn^/mf; With 

^' = v^/vl, model /, 

= —1, model II. 
When C = C = 0, the above result (H) reduces to that of SM case||20|, |21 . 



The renormahzation group equation satisfied by the coefficient functions Ci(/i) is 

/^|;^^(/^) = E(7n.A(/^)- (7) 

Where the anomalous dimension matrix is calculated in practice by requiring renormahzation 
group equations for Green functions with insertions of composite operators to be satisfied order 
by order in perturbation theory. 

After evaluating the loop diagrams, we find the following leading order weak mixing of 
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operators, with the Q, R part agrees with ref . ||20[| . 
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These mixing are all between operators induced by tree-diagrams and operators induced by 
loop-diagrams. The vanishing log(yU^/m^) terms in previous matching are regenerated here by 
renormalization group equation (|^. 
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The QCD anomalous dimensions for each of the operators in our basis are 
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The mixing elements between the two matrix (P) and ([1^) are all zero in leading log approx- 
imation. There are some differences in the anomalous dimension matrix comparing with Cho 



and Grinstein's result pO[. If we omit a symmetric factor of 1/2 in calculating Feynman diagram 



like Fig.3, the results will agree with them|21]. But according to similar diagram calculation like 



ref.[^, |T0|, this factor can not be omitted. For instance, diagram like Fig. 3 is also appeared in 
calculation of anomalous dimension of the gluon magnetic moment-type operator Og in ref.[^, |10| 
with a symmetric factor 1/2. This is also shown in general Feynman gauge calculation, and it 
give results same as the above. After these changes, the whole matrix can be easily diagonalized, 
and gives all real eigenvalues, which is required to maintain hermiticity of the effective hamilto- 
nian at all renormalization scales. While in ref.||20|, it can not. In their case, some eigenvalues 
are complex. 

The solution to renormalization group equation appears in obvious matrix notation as 



C(/i2 



exp / dg—— 

53(m) f3{g) 



'111 



After inserting anomalous dimension (|8|-|T0D, we can have the coefficients of operators at /i = M^. 



3 QCD Corrections from fi = Mw to = mh Scale 

In order to continue running the basis operator coefficients down to lower scales, one must 
integrate out the weak gauge bosons and would-be Goldstone bosons at /i = My/ scale. The 
diagrams are displayed in Fig.4. From the second and third matching equations of this figure, one 
finds the following relations between coefficient functions just below(-) and above(+) yU = Mw'- 

Coi^iM^) = C^oi,(M+), 
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In addition to these, there are new four-quark operators from the first equation of Fig.4p, H, |T0 
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with coefficients 



Ci{Mw) = 0, ^ = 1,3,4,5,6, C2{Mw) = l. 



To simphfy the calculation and compare with the previous results, equations of motion(EOM) [ 
is used to reduce all the remaining two-quark operators to the gluon and photon magnetic mo- 
ll 



ment operators 0|^^ and 0|j:j. To be comparable with previous results without QCD corrections 
from nitop to M]y, operators O^^j^, are rewritten as O7, Og hke ref. p!0[| , 
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The operator basis now consists of 8 operators. The effective Hamiltonian appears just 
below the W-scale as 
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For completeness, the explicit expressions of the coefficient of operator Os and O7 at /i = 
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Since they are expressed by coefficients of operators aX ji = rrit and QCD coupling a^, it is 
convenient to utilize these formula. 

If the QCD corrections from rritop to Mw are ignored (by setting asirrit) = as{Mw) in 
eqn.([T7|),([l8|) ), the above results(p^ (p^Sf) would reduce to the previous results^, |15| exactly, 
where the top quark and W bosons are integrated out together: 

CriMw) = -l-Ai^)-l(^-^yA{y)+eB{y) (19) 



Cs{Mw) = -l-D{x)-l(^yD{y)+eE{y), (20) 



2 ' ' 6 \vi. 

2 ^ ^ 6 \vi, 
with A(x), B(y), D(x), E(y) defined in ref.[§. 

The effects of QCD corrections to Cj{My[r) and Cs{Myy) can easily be seen from Fig. 5 and 
Fig. 6. Here W boson mass is taken as Mw = 80.22GeV, The top quark mass rrit = 174GeV, and 
the QCD scale is taken as A{^^Jd = 175 MeV^, Mh+- = SOOGeV. The resuhs of model I is 
displayed in Fig. 5. Except for small values of V2/V1, Ci{M]y) and Cs{M]y) are both suppressed 
by QCD corrections from rrit to Mw- At V2/V1 = 10, Ci{Mw) is suppressed 17% and Cs{Mw) 
12%. Fig. 6 gives results calculated in model II, the differences are that |C7(Mvk)| and |C8(Mvk)| 
are both enhanced by QCD corrections for all values of V2/V1. At small values of ^2/^1, |C7(Mvi/)| 
is enhanced 10%, |C8(Mvi/)| is enhanced 8%. Since C-j{Mw) and Cii{Mw) are both the input of 
the following QCD running from Mw to mf,. It is expected to change the final result. 

The running of the coefficients of operators from ^ = Mw to fi = rrih was well described 



in ref. g, [T0[- After this running we have the coefficients of operators at = nib scale. Here 
rub = 4.9GeV is used. Except small values of V2/V1, both Cjinib) and C^ijnb) calculated in 
model I are suppressed in comparison to values obtained by ref.p, 53, where the QCD running 



from rrit to Mw is neglected. While in model II both \C-j{mb)\ and \Cs{Tnb)\ are enhanced. 
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4 The B Xgj decay rate 



The leading order b —>■ sj matrix element of -f^e// is given by the sum of operators O5, Oq and 



O7 in our effective theoryP, 10 



<H,ff>=-2V2GFV*Vt,{CM + Qd[C,{fi)+3Ce{fi)]}<\07\>. (21) 

Therefore, the sought amplitude will be proportional to the squared modulus of 

C'/^im,) = Cr{m,) + [C,{m,) + 3CQ{m,)] (22) 

instead of |C7(m;,)p itself. 
Following ref.[|, |, 0, 

BR(B Xs-i)/BR{B XcCV) ~ T{b s-f)/T{b ceV). (23) 



Then applying eqs.(pTD,(P2D, one finds 



-1 



. (24) 

Bn{B — > XcCv) -ngimc/mi,) \ in j 

where g{mc/mh) — 0.45 and flmc/mh) ~ 2.4 corresponding to the phase space factor and the 
one-loop QCD correction to the semileptonic decay, respectively [^. The electromagnetic fine 
structure constant evaluated at the b quark scale takes value as aQED{f^b) = 1/132.7. Afterwards 
one obtains the B decay rate normalized to the quite well established semileptonic 

decay rate. If we take experimental result Br{B XcCV) = 10.8%|2^, the branching ratios of 
B — > Xs'j is found. 

The decay results are summarized in Fig. 7 and Fig. 8 as functions of V2/V1, with different 
charged Higgs mass 170GeV, 300GeV, 600GeV and 900GeV. The CLEO upper limit is also 
shown solid line. 
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For 2-Higgs doublet model I, the decay rates including QCD running from rrit to Mw are 
suppressed except small values of V2/V1. At V2/V1 = 10, 171^+- = SOOGeV, the suppression is 
23%. As V2/V1 — > 0, this kind of model goes back to SM. In this small V2/V1 region, the QCD 



corrections give 7% enhancement which corresponding to SM case||2T|. Therefore, the restrictions 
to this model parameters are less tight than previous predictions especially for large values of 
V2/V1. In Fig. 7, one can see that, for lower mass of charged Higgs, small values of V2/V1 are still 
allowed. The parameter space is still open. The decay channel t ^ hip can still exist for a wide 
region. 

The decay rates in 2-Higgs doublet model II are enhanced up to 5% more than previous 
calculations when charged Higgs mass is lower than 800 GeV. Although this percentage is not 
very large, the absolute values of the decay rates are sure higher. In Fig.8, one can easily see 
that, as charged Higgs mass lower than 700GeV, all values of V2/V1 are excluded by CLEO's 
b ^ experiments. For larger charged Higgs mass, small value of V2/V1 is still allowed. It is 
obvious that, the top decay channel t ^ bcj) is already ruled out in this kind of 2-Higgs doublet 
model. 

In the Supersymmetry model, the Higgs sector is the same as model II, but the large chargino 



contribution cancels much the charged Higgs contribution!^. There is sure a suppression of 
decay rates from 2-Higgs doublet model II. The magnitude is more complicated, since it depends 
on various Supersymmetry models and chargino mass. When chargino mass is of order 0{mtop) 
or higher, it should be integrated out at /i = rritop] if it is of order 0{Mw) or lower, it should 
be integrated out at /i = Mw scale. Further more, for a complete Supersymmetry calculation. 



there are also contributions from other SuSy particles which need a detailed calculation p7| 
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5 Conclusion 



As a conclusion, we have given the full leading log QCD corrections (including QCD running from 
rritop to M\y) to 6 — s7 decay in 2-kind 2-Higgs doublet models. 

The QCD running from to M^r suppresses the b ^ s'y decay rate in model I, and enhances 
the decay rate in model II. 

Although this result is not quite different from the previous calculations, our improvements 
lie in reducing some theoretical uncertainties. After these changes, restrictions from 6 — > 57 
decay to 2-IIiggs doublet model parameters are less strict in model I and more tight in model II 
than previous predictions. It is shown that the decay 6 — > 57 is by far the most restrictive process 
in constraining the parameters of the charged Higgs boson sector in 2-Higgs doublet model. 

Note added: 

After this paper was submitted for publication, the paper|2^ by H. Anlauf came to our 



attention; the author studied the supersymmetric contributions as well as charged Higgs contri- 
butions in 2HD model II. Our results consist of 2-kinds of 2HD models; and agree with their 
result of model II at leading order (Mpi//mj, mH/rrit). Further more, we also included high order 



contribution which can not be neglected(In ref. this was picked up later after running to M] 



iw 



scale to match onto results without QCD running from rrit to Mw)- 
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Figure Captions 

Fig.l Leading order matching conditions at the top quark scale for the IPI Green functions in 
the full theory and in the intermediate effective field theory. 

Fig. 2 One loop matching conditions at the top quark scale for the IPI Green functions in the 
full theory and in the intermediate effective field theory. 

Fig. 3 One of the Feynman diagram in calculating anomalous dimensions, with the heavy dot 
denoting high dimension operator. 

Fig. 4 Matching conditions at = M]y for four quarks and two quarks IPI Green functions in 
the intermediate effective field theory and effective field theory below W scale. 

Fig. 5 The photon and gluon magnetic moment operator's coefficients C7{Mw) and Cs{Mw) in 
model I. The solid lines are our QCD corrected results, the other two are uncorrected ones. 

Fig. 6 The photon and gluon magnetic moment operator's coefficients Ct{Mw) and Cs{Mw) in 
model II. The solid hues are our QCD corrected results, the other two are uncorrected ones. 
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Fig. 7 BR{B Xgj) of model I as function of V2/V1 for different charged Higgs masses. The 
sohd line is the upper limit of CLEO. This line upper is the excluded region. 

Fig. 8 BR{B — > Xg^) in model II as function of V2/V1 for different charged Higgs masses. The 
solid hne is the upper hmit of CLEO. 
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